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Abstract 

We made high-resolution spectroscopic observations of hmb-spicules in Ha using the Vertical 
Spectrograph of Domeless Solar Telescope at Hida Observatory. While more than half of the observed 
spicules have Gaussian line-profiles, some spicules have distinctly asymmetric profiles which can be fitted 
with two Gaussian components. The faster of these components has radial velocities of 10-40 km 
and Doppler-widths of ^0.4 A which suggest that it is from a single spicule oriented nearly along the 
line-of-sight. Profiles of the slower components and the single- Gaussian type show very similar character- 
istics. Their radial velocities are less than 10 km s~^ and the Doppler-widths are 0.6-0.9 A. Non-thermal 
"macroturbulent" velocities of order 30 km s~^ are required to explain these width-values. 

Key words: line:profiles — Sun:chromosphere — Sun:spicules 



> 

in 
\o 
m 

(N 
O 
O 



1. Introduction 

Spicules are one of the fundamental elements of the 
quiet solar chromosphere, so it is important to know the 
physical properties of spicules to understand the quiet 
Sun. Pasachoff, Noyes, Beckers (1968) did a pioneer- 
ing observation of limb spicules with many chromospheric 
emission lines, and derived the basic morphological and 
spectroscopic properties. They also suggested the rota- 
tion of spicule gas from the inclined spectra. Krat & 
Krat (1971) also made spectroscopic observations of limb 
spicules in 5 spectral lines, and reported that emission- 
line profiles are classified into narrow ones and wide ones. 
However, they were not sure whether they observed sin- 
gle individual spicules, since the spatial resolution of their 
study was several arcseconds. Summarizing these obser- 
vational works, Beckers (1972) reviewed the morphology 
and spectroscopic properties of spicules and summarized 
that a spicule has a diameter of 400-1500 km and mean 
Doppler velocity of 10 km s~^. Kulidzanishvili & Nikolsky 
(1978) observed 650 Ha line profiles of 25 spicules, and 
categorized them into two groups. One group is of rel- 
atively small intensity and narrow emission profiles, and 
the other is of brighter and wider ones. The character- 
istics of the latter group was thought to be due to the 
overlapping of unresolved spicules. They suggested that 
the Ha spicules have temperature of 6000 K and non- 
thermal velocities of ^25 km s^^. Follow-up spectro- 
scopic observations of limb spicules under much better 



seeing conditions were done by Kulidzanishvili (1980) with 
the 53-cm Lyot coronagraph of Abastumani Astrophysical 
Observatory and by Hasan & Keil (1984) with Sacramento 
Peak VTT. From the study of line-of-sight (LOS) veloc- 
ity and line-width distribution along the axis of spicules, 
they got a result that, in majority of spicules, the Ha line- 
width does not change spatially along the axes of spicules. 
Kulidzanishvili (1980) suggested that the line width of the 
majority of spicules are due to the overlapping of gas mov- 
ing with relative velocities of 20-30 km s~^. 

Nishikawa (1988) used Ha filtergrams to determine the 
diameter and the motion of limb spicules. This morpho- 
logical study found that the speed of the apparent up-and- 
down motion of large spicules is ^50 km s~^. Assuming 
ballistic motion, he estimated the initial velocity to be 
around 100 km s^^, though there has been no other report 
of such high velocities. The typical diameter of spicules 
was about 500 km, although the best images show thinner 
components of 200 and 350 km in diameter. 

Recent studies using Call H filtergram from Hinode 
has shown that spicules move more dynamically and have 
much smaller diameters than previously thought. De 
Pontieu et al. (2007b) found that the apparent maximum 
velocity of "type I" spicules is 15-35 km s~^, and the ap- 
parent upward velocity of "Type H" spicules is 40-300 
km s~^. Suematsu et al. (2008) reported that a spicule 
consists of highly dynamic multi-threads as thin as a few 
tenths of an arcseconds and shows lateral movement or 
oscillation with rotation. 
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Although Hinode's direct imaging provides us highly 
valuable images, we need spectroscopic observations of 
spicules with high spatial and temporal resolutions to 
know the actual physical state of the spicular gas. As 
Sterling (1998) pointed out, knowledge of the actual phys- 
ical properties is essential to theoretically clarify the 
dynamics and ejection process of spicules with numeri- 
cal simulations. Recently, Pasachoff, Jacobson, Sterling 
(2009) presented the ground-based observations with very 
high spatial resolution done at the Swedish 1-m Solar 
Telescope on La Palma. They analyzed the imaging ob- 
servations at five wavelengths around Ha. Although their 
observation provides us a new statistics on morphological 
and dynamical properties of limb spicules, they have de- 
rived the spectroscopic quantities such as LOS velocities 
under the rather simplified assumption of the Gaussian 
emission profiles. At present, further analyses of full line 
profiles in detail remain to be done to find the thermody- 
namical properties of spicular gas. 

In our series of papers, we present the analyses of spec- 
troscopic observations of limb spicules in Ha with high 
spatial resolution, which is comparable to the best reso- 
lution of ground-based observations. This paper reports 
and discusses the line widths and the radial velocities of 
spicules derived from our spectral images. 

2. Observation 

The observations were made on September 8, 2005 at 
Hida Observatory, using the Vertical Spectrograph of the 
Domeless Solar Telescope. The 3rd order spectra gave 
us the dispersion of 0.300 A mm^^ and the spectral res- 
olution 0.0123 A (which corresponds 0.56 km s~^) with 
the slit width 0.1 mm. Ha spectra were taken by Kodak 
Megaplus 1.6i CCD camera with 150 msec exposure time. 
The slit position was fixed during each series of observa- 
tions and the camera continuously took the spectral im- 
ages with a 1 sec time cadence. 

We targeted the north polar limb which was covered 
by a large coronal hole at that time (figure 1). The slit 
was positioned parallel to the limb and at the tops of the 
spicules to get the optically thin profiles (figure 2). We 
were able to obtain more than 3500 frames of Ha spectral 
images with high spatial resolution during two series of 
observations around 01 UT and 06 UT. The height of the 
slit position at the nearest point from the solar limb seen 
in Ha-hO.9 A was 4200 km for 01 UT and 3700 km for 06 
UT. 

Among the numerous spectral images, we first discarded 
the frames showing dark lane of self-absorption around the 
Ha line center, since their profiles with fiat-bottomed dip 
require the radiative transfer analysis to interpret and we 
wish to concentrate our profile analysis to optically thin 
emissions. Then, we further searched the frames showing 
good spatial resolution. As a result, seven frames shown 
in figure 3 were picked up for analysis, three from the 
sequence around 01 UT and four from 06 UT. The spatial 
FWHM values of the sharp streaks measured at ±0.8 A 
are 1000-1900 km. 




□T 195 20050902-191347 



Fig. 1. SOHO EIT shows a coronal hole around the north 
pole. The rectangle indicates the observed region. 
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Fig. 2. Ha images of the observed regions at 01 UT and 06 
UT. The horizontal line is the spectrograph slit. Two vertical 
lines are fiducial hairlines. The images are presented in loga- 
rithmic intensity scale to enhance the details around the slit. 
The arrows indicate the positions of two exceptional events 
mentioned in section 5. 
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3. Ha spectra of spicules analyzed in this paper. Two vertical lines correspond to the fiducial hairlines. 
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Recent Hinode/SOT observations found that the spa- 
tial widths of spicules are around 200 km or less. The 
spatial resolution of our spectroscopic observation is much 
coarser than the true widths of spicules. Unfortunately, 
our observation was before the launch of Hinode and we 
could not make use of the high-resolution images of SOT 
to compare. 

3. Data Analysis 

Dark subtraction and flat fielding were applied to the 
observed spectral images, then a mean sky spectrum was 
subtracted to obtain the proper emission profile of the 
spicules. 

3.1. Method of Profile Fitting 

We got the profiles at every pixel position along the 
slit between the hairlines in figure 2. The spacing is 0.15 
arcscc per pixel, so we got about 650 profiles from each 
frame. Telluric lines were removed from the profiles, then 
each profile was fitted first with a single Gaussian curve as 
a function of wavelength from Ha line-center (AA), which 
determines relative peak intensity /q, Doppler width wy>, 
and Doppler shift AAd of the emission. 

/,ingic(AA)=/oe \^^) (1) 

Some profiles were asymmetrical and were not well fit- 
ted with a single Gaussian curve. The goodness of fit was 
eveluated by a reduced chi-square (x^) statistic. First 
we eliminated the low-intensity (/q < 3.0 sky brightness) 
profiles as they had a low signal-to-noise ratio. Then we 
introduced xl ^ threshold for a good fit, because the 
effect of the values are influenced by intensity. For the 
profiles with a large x/Io value, we tried to fit the ob- 
served profile with a sum of two Gaussian curves, and get 
two sets of parameters. 

/double (A A) /slow e ^ "^low ; +/faste V "fast ; (2) 

The subscripts "slow" and "fast" denote the component 
with smaller and larger Doppler-shift values respectively. 

Figure 4 is an example of a single-peaked and almost 
symmetrical profile which requires the double Gaussian 
fitting. While there remains substantial deviations at the 
peak and far wing of the profile in the case of the single 
Gaussian fitting (left panel), the double Gaussian fitting 
(right panel) can well follow the whole curve of the ob- 
served profile. The goodness of fit (chisq) becomes 10 
times better in the case of the double Gaussian fitting. 
Similarly, for single-peaked but fairly asymmetrical pro- 
files, the double Gaussian fitting works very well in our 
data. 

All of the analyzed profiles are successfully fitted with 
single or double Gaussian curves. This confirms that the 
emissions in our selected frames were from optically thin 
part of the spicules. 



3.2. Effect of Atmospheric Image Blurring 

The spatial pattern of the spicule emission is blurred 
by the seeing effect. We tried to recover un-blurred spec- 
tral images with the help of 1-D Wiener filter method 
and checked the effect of atmospheric image blurring on 
the spectral parameters of emission profiles. We assumed 
that the point-spread-function (PSF) of the atmospheric 
blurring is a normal distribution. The width of the PSF 
was estimated from the smallest spatial width of observed 
emission and assumed to be uniform over the field of 
view. The power spectrum of the noise was estimated from 
the observed intensity fluctuation at the continuum wave- 
length domain where no spicular emissions were detected. 
Since its power spectrum was nearly flat, we assumed the 
noise to be white-noise. 

An example of compensated spectral image is shown in 
figure 5 with its original observed image. We can see that 
the emission streaks became sharper in spatial direction 
and more clear in the "deconvolved" image. 

The same spectral fitting method as described in the 
subsection 3.1 was applied to the deconvolved image to see 
the effect of atmospheric blurring on spectroscopic quan- 
tities. Figure 6 shows the results for two typical profiles. 
We can see that the spectral parameters do not change 
so much between the originally observed and the decon- 
volved image, except the intensity of the fast component. 
We performed the same comparison on all the profiles in 
the deconvolved spectral image, and estimated the effect 
of the atmospheric blurring on the spectroscopic parame- 
ters. For the single-Gaussian fitted profiles, /q changes less 
than ±5%, and ujd changes less than ±1%. The change 
of AAd is smaller then 0.05 A. In the case of asymmet- 
ric profiles, the changes of lUsiow and AAgiow are also very 
small. Wfast changes ±8%, and AAfast changes ±6%. /siow 
decreases 4.5%, with 9.6% deviation. Only /fast shows 
significant change, with average increase of about 14%. 
Increase in /fast is expected since the deconvolution con- 
verges the blurring of isolated emission. Thus the spec- 
troscopic parameters, especially the widths and the shifts 
of fitted components are not so influenced by the atmo- 
spheric blurring of the images. 

Our method of seeing correction is based on a few as- 
sumptions, which are (1) the uniformity of seeing over the 
field of view, (2) whiteness of the power of background 
noise, etc. In the deconvolved image, there appear faint 
and equally-spaced emission streaks in the featureless part 
of the original spectral image. These equally-spaced emis- 
sion patterns are probably artificial due to the drawback 
of the method we used. So we will proceed our discus- 
sion based on the analysis without the correction of at- 
mospheric blurring. 

4. Results 

Line profiles of spicules can be grouped into three types. 
Profiles of the first type (Type A) can be fitted well with 
a single Gaussian curve. The second type (Type B) is 
not well fitted with single Gaussian, but it is nearly sym- 
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Fig. 4. Comparison of results with single and double Gaussian fitting. Though the original emission profile is almost symmetrical, 
single Gaussian fitting (left) cannot reproduce the observed profile. Double Gaussian fitting for the same profile (right) works well. 
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Fig. 5. The observed spectral image (left) and the processed image with the Wiener deconvolution method (right). 



M. Shoji et al. 



[Vol. 



c 

01 



m 

> 



01 



original 




'.2^i509D8J3 10608 position=573 

deconvolved 



Single 

chisq=0.005S 



-3 -2 -1 J 
wavelength (A) 



I: 

iX: 



5.8182 

0.0951 
0.7027 




-2-1 \ 
wavelength (A) 



Single 

chisq=0.00&6 



I: 

AA: 
w: 



5.8266 

0.1028 
0.6908 



m 
r. 



01 




200509 DS_0 10608 



Double 

chisq-0.0071 

Sfow 

4.6S25 
0.0734 
0.6920 

1 .9809 
-0.6666 
0.4262 



-J5. -2-10 1 
wavelength (A) 



position =552 
deconvolved 



c 




-3-2-10 1 
wavelength (A) 



Double 

chisq=0,0077 

Slow 

I: 4.4248 

AA: 0.0945 

w: 0.6727 



Fast 

I: 

w: 



2.4219 
-0-6542 
0.4552 



Fig. 6. Examples of the emission profiles to compare between the original and the "deconvolved" spectral images shown in figure 5. 
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Fig. 7. Examples of the line profiles. The unit for wavelength shift (AA) and Doppler width values is A. (a) Typical Type A 
profiles, which are well fitted with a single Gaussian curve, (b) Typical Type B profiles, which are symmetrical, but cannot be fitted 
well with a single Gaussian curve. 
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Fig. 8. Examples of the typical line profiles of Type C, which are asymmetrical and separated to slow and fast Gaussian components. 
The unit is the same as figure 7. 
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Fig. 9. (a) An example image sliowing Dopplor shifted streaks of spicules, (b) Radial velocities determined by single or double 
Gaussian fitting are plotted over the same image as (a). Error bar (light blue) indicates the standard deviation of each velocity 
value, (c) Doppler widths of the profiles at the same positions as (b). Note that two values are plotted at the same axis position 
when the profile is fitted with double Gaussian curves. 
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Table 1. Classification of Ha profiles. Number of points 
where the line profile is categorized into three types are tabu- 
lated. Type A: well fitted with single Gaussian profile. Type 
B: ill fitted with single Gaussian, but almost symmetrical. 
Type C: asymmetrical profile which can be separated into 
slow and fast components. 



Time (UT) 


Type A 


TypeB 


TypeC 


01:06:08 


226 


96 


48 


01:11:12 


234 


149 


75 


01:24:22 


397 


78 


36 


06:10:45 


369 


161 


53 


06:12:50 


351 


111 


69 


06:16:44 


253 


103 


137 


06:23:34 


233 


128 


63 


Total 


2063 (61%) 


826 (25%) 


481 (14%) 



metrical around the Ha line center (figure 4). The last 
type (Type C) is asymmetrical, and separated into an 
almost stationary component and a highly shifted compo- 
nent. Note that the difference between Types B and C 
is actually a matter of degree, and there are many pro- 
files in-between. For argument's sake, we take the profiles 
with I AAfast/AAsiowl > 4.0 as Type C and treat all other 
double Gaussian profiles as Type B. Figures 7 and 8 show 
the examples of typical profiles of each type, and table 1 
lists the number of points categorized into each type. 

Figure 9 is an example showing estimated LOS veloci- 
ties and widths with the original spectral image. In figure 
10, we summerize the distribution of Doppler widths, ra- 
dial velocities, and equivalent widths in histogram format. 

4.1. Doppler Width 

Figure 10a shows that the Doppler widths of the Type A 
profile (wd) are 0.6-0.9 A. Widths of the slow components 
of the Type C profile {wgiow) are also 0.6-0.9 A, while 
widths of the fast component (wfast) are 0.3-0.6 A. 

To identify the line-broadening mechanism, we calcu- 
lated the theoretical Doppler width of Ha profile with 
given temperature T, microturbulence ^, and macrotur- 
bulence S, 



A 2kT 



w -■ 



-^2 + 22 



(3) 



where A is the wavelength of Ha line, c is the light velocity, 
k is Boltzmann constant, and m is hydrogen mass. 

The result is shown in figure 11. When S = 0, the ob- 
served value Wfast = 0.35 A is found within a reasonable 
range of temperature and microturbulence, T ~ 15000 K 
and ^ ~ 5 km s~^. On the other hand, even if we assume 
as high a temperature as T = 20000 K with microturbu- 
lence ^ = 20 km s~^ we cannot get a width more than 0.6 
A. While Beckers (1972) listed the most likely tempera- 
ture as 9000-16000 K, Makita (2003) derived ionization 
temperature for hydrogen to be 5200 K and excitation 
temperature for Ha to be 5020 K, and other authors also 
reported lower temperatures. Therefore, it is unlikely that 
the tempetrature of spicules to be higher than 20000 K. 
And considering that microturbulent velocity cannot be 
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Fig. 10. Histograms made from the results of all frames. 
"Slow" and "fast" mean the slow and the fast components 
of the Type C profile, (a) Doppler width, (b) Radial velocity 
(absolute value), (c) Equivalent width. 
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Fig. 11. Theoretical Dopplcr width, calculated as a function 
of temperature with different values of microturbulent veloc- 
ity { (km s~^) and macroturbulent velocity. 

supersonic, we need additional "macroturbulence" to ex- 
plain the observed Wsiow and w\y . 

From figure 11, we find that the required S for Wsiow and 
uiD is about 30 km s~^. Here we use the term "macro- 
turbulence" to include any non-thermal random velocities 
other than microscopic turbulence. There exists many 
spicules with various LOS velocities along the LOS of ob- 
servation. Superposition of their emissions will result in 
the broadening of the observed profile. "Macroturbulent 
velocity" is a measure of the dispersion of their LOS ve- 
locities. 

J^.2. Radial Velocity 

Figure 10b is the histogram of radial (line-of sight) ve- 
locities. With our definition for Type C, the distribution 
of the fast component's parameters hardly overlaps with 
that of the slow component's. Radial velocities derived 
from Type A profiles are less than 10 km s~^, and the pro- 
files of the slow component of Type C give similar values. 
The distribution of the fast component peaks at around 
25 km s"-'^ and there is also a small peak at around 40 
km s^^. 

4-3. Correlation between Width and Velocity 

Figure 12 is a correlation plot between radial velocities 
and Doppler widths for Types A and C. Again, it is clear 
that the slow component of the Type C profile has a simi- 
lar distribution to the Type A profile. The velocity of the 
fast component spreads from 10 km s~^ to 40 km s~^, but 
its width is independent of radial velocity. 




10 20 30 40 50 
velocity (km/s) 



Type C 




10 20 SO 40 50 
velocity (km/s) 



Fig. 12. Correlation between velocity (absolute value) and 
Doppler width. Crosses are for the slower component and 
diamonds are for the faster component. Different color repre- 
sents different time. 

5. Summary and Discussion 

We made high-resolution spectroscopic observations of 
limb spicules in Ha. While more than half of the ob- 
served spicules have Gaussian line profiles (Type A) , some 
spicules have distinctly asymmetric profiles which can be 
fitted with two Gaussian components (Type C). The faster 
component of Type C profiles have radial velocities of 10- 
40 km and Doppler width of ^0.4 A. This width can 
be explained with temperature ~ 15000 K and microtur- 
bulence ~5 km s^^. On the other hand, widths of Type A 
profiles and the slower component of Type C profiles are 
0.6-0.9 A. Macroturbulent velocities of order 30 km s~^ 
are required to explain this large width. 

Now we will discuss the physical interpretations of these 
results. 

Figure 2 shows that the slit-line intersects many 
spicules, and it is natural that the LOS intersects sev- 
eral spicules, too. Though it is difficult to observe a sin- 
gle spicule without superposition of other spicules. Type 
C profile enables us to separate the emission of single 
spicule as the fast component of double-Gaussian fitted 
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components. Line width of the fast component can be ex- 
plained with temperature ~15000 K and microturbulence 
~5 km s~^. Broad width of the Type A profiles suggest 
that many spicules along the LOS contribute to the ob- 
served emission. Figure 10c shows that their equivalent 
widths are five times or more larger than those of the fast 
component. 

The LOS velocity 20-40 km s~^ of the fast compo- 
nent is consistent with the velocity of the apparent up- 
and-down motion of spicules reported by many authors 
(e.g. Nishikawa 1988; Dc Pontieu et al. 2007b; Pasachoff, 
Jacobson, Sterling 2009). Because most spicules are 
nearly vertical to the solar limb (Beckers 1972), their LOS 
velocities are much smaller than their actual velocities. 
We suppose that the profiles of Type A and slow compo- 
nent of Type C are from these nearly vertical spicules, and 
that the fast component is from the few spicules which are 
almost parallel to the LOS. 

As shown in section 4.1, we need additional macro- 
turbulent i.e. non-thermal random velocities of order 30 
km s~^ to explain the Doppler width of the slow compo- 
nent of Type C and the Type A profiles. Other obser- 
vations also reported the high non-thermal broadening. 
For example, Makita (2003) suggested that spicules have 
a turbulence of ~20 km in the active region, based 
on the analysis of Ca II H and K profiles in the flash spec- 
trum of the 1958 eclipse. And Mariska, Feldman, Doschek 
(1978) reported average non-thermal velocity of 28 km s~-^ 
in EUV line profiles at the quiet limb. Next we proceed 
to discuss the source of macroturbulence. 

One of the possible sources is ejecting motion of spicules 
inclined toward or away from the observer. Variations 
in speed and inclination of ejection produce the disper- 
sion of LOS velocities. As discussed above, most spicules 
have small inclinations and so their LOS velocities will be 
much smaller than 10-40 km s~^. Distributions of appar- 
ent inclinations reported by Pasachoff, Jacobson, Sterling 
(2009) have peaks at 10 and 25 degrees, in which case 
the LOS velocity of ejection speed 40 km s~^ will be at 
most 7 and 17 km s^^. They are too small to explain 
the observed widths. There is also the height difference of 
ejection speed. The spicules nearer to us or farther than 
the limb of the sun would be observed in superposition 
against lower heights of spicules whose bases are on the 
limb. However, the speed at the tops is no faster than the 
lower part (Sterling 1998; De Pontieu et al. 2007b), and 
it will contribute little to the line broadening, similar to 
the argument above. 

Another possible source is the lateral motion of spicules 
due to the Alfven wave disturbance recently found in 
Hinode/SOT Call H fihergrams (De Pontieu et al. 2007a; 
Suematsu et al. 2008; He et al. 2009b), or the kink wave 
observed by He et al. (2009a). De Pontieu et al. (2007a) 
showed that the distribution of transverse displacements 
of the spicules agrees with the velocity amplitudes around 
20 km s""'^, while He et al. (2009b) reported the velocity 
amplitude of high-frequency Alfven waves to be 4.7-20.8 
km s"'^. Velocity amplitude of the kink wave reported by 
He et al. (2009a) is less than 8 km s^^. Superposition 



effects of these motions can broaden the line profiles. 
Because spicules are believed to be ejected along magnetic 
field lines, transverse Alfven waves broaden the Ha line 
profiles of Type A or slow component of Type C. MHD 
simulation also showed that the amplitude of Alfven waves 
at the chromospheric height to be ~20 km s~^ (Suzuki & 
Inutsuka 2006). However, when the sinusoidal motions of 
independently disturbed spicules are superposed, its stan- 
dard deviation will be 1/ \/2 of amplitude. Thus oscillation 
of velocity amplitude 20 km s~^ will contribute only 14 
km s~^ to the macroturbulcnt velocity. As for the fast 
components, we assume that their surrounding magnetic 
field is nearly parallel to the LOS, so the Alfven waves 
will not contribute to line-broadening, agreeing with their 
narrow widths even if they are composed of multiple finer 
threads. Kulidzanishvili (1980) reported the oscillations 
of the radial velocities of limb spicules with the period 
of 3-7 min. However, their amplitudes were less than 10 
km s~^. which arc too small to contribute to our observed 
line- widths. 

Even when we include superposition effects of ejecting 
motions and Alfven wave broadening, it is still insufficient 
to explain the macroturbulcnt velocities of ~30 km s^^. 
So there must be some additional but unidentified sources 
for the non-thermal broadening of spicule emission pro- 
files. Alfven waves with higher frequencies and shorter 
wavelengths than those observed until now may be present 
and contribute to the broadening. 

Let us mention two interesting events found in our anal- 
ysis. Their characteristics are different from the common 
spicules discussed so far. While most of radial velocities 
are less than 40 km s~^, these two cases have much larger 
( > 45 km s~^) radial velocities as seen in the upper right 
area of figure 12. Their Doppler widths are exception- 
ally large, too. They are in the spectra at 06:16:44 UT 
and 06:23:34 UT, and their faint and broad profiles can be 
found in figure 8. The positions of these events are marked 
in figure 2 where there is a gap in the bush of spicules, and 
where the slit is the nearest to the limb. We need differ- 
ent broadening mechanism for these events. It may be an 
event with a magnetic reconnection, or something related 
to the spicule formation. This phenomenon may corre- 
spond to the Type H spicules of De Pontieu et al. (2007b) 
or to macrospicules (Bohlin et al. 1975; Yamauchi et al. 
2003) or to polar-jets (Shibata et al. 1992; Shimojo et al. 
1996; Savcheva et al. 2007). The relations are yet to be 
ascertained in future studies. 

In this paper, we treated the results statistically with- 
out regard to their positions. However, figure 3 shows that 
some spicules appear as tilted streaks. This is possibly an 
evidence of rotation of spicules suggested by Pasachoff, 
Noyes, Beckers (1968), and its spatially blurred emission 
may contribulte to broaden the line profiles. 

Some physical values are not determined by the Ha pro- 
file alone. In order to study the true origin of macrotur- 
bulence or of the high-velocity events, we need simulta- 
neous spectroscopic observations of spicules in multiple 
spectral lines, which can be achieved with the Horizontal 
Spectrograph at Hida Observatory. We also need high 
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spatial-resolution images of Hinode/SOT simultaneously 
obtained with the spectroscopic observation to understand 
the details of spicules. We will plan our next observation 
to accommodate to those requirements. 
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